WEAKLY DIRECTED SELF- AVOIDING WALKS 



AXEL BACHER AND MIREILLE BOUSQUET-MELOU 

Abstract. We define a new family of self-avoiding walks (SAW) on the square lattice, called 
weakly directed walks. These walks have a simple characterization in terms of the irreducible 
bridges that compose them. We determine their generating function. This series has a complex 
singularity structure and in particular, is not D-finite. The growth constant is approximately 
2.54 and is thus larger than that of all natural families of SAW enumerated so far (but smaller 
than that of general SAW, which is about 2.64). We also prove that the end-to-end distance 
of weakly directed walks grows linearly. Finally, we study a diagonal variant of this model. 



1. Introduction 

A lattice walk is self-avoiding if it never visits the same vertex twice (Fig. [T). Self-avoiding 
walks (SAW) have attracted interest for decades, first in statistical physics, where they are 
considered as polymer models, and then in combinatorics and in probability theory [55]. However, 
their properties remain poorly understood in low dimension, despite the existence of remarkable 
conjectures. See |25| for dimension 5 and above, and [7] for recent progresses in 4 dimensions. 




Figure 1. A self-avoiding walk on the square lattice, and a random SAW of 
length 1,000,000, constructed by Kennedy using a pivot algorithm |22| . 

On two-dimensional lattices, it is strongly believed that the number c„ of n-step SAW and 
the average end-to-end distance D n of these walks satisfy 

c n ~ a/i n n 7 and D n ~ nn v (1) 

where 7 = 11/32 and v = 3/4. Several independent, but so far not completely rigorous methods 
predict these values, like numerical studies HH], comparisons with other models [5J 125] . 
probabilistic arguments involving SLE processes [24], enumeration of SAW on random planar 
lattices |13|... The growth constant (or connective constant) fx is lattice-dependent. It has recently 

been proved to be v2 + \/2 for the honeycomb lattice [T2] , as predicted for almost 30 years, and 
might be another bi-quadratic number (approximately 2.64) for the square lattice |21| . 

Given the difficulty of the problem, the study of restricted classes of SAW is natural, and 
probably as old as the interest in SAW itself. The rule of this game is to design new classes of 
SAW that have both: 
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- a natural description (to be conceptually pleasant), 

- some structure (so that the walks can be counted, and their asymptotic properties de- 



This paper fits in this program: we define and count a new large class of SAW, called weakly 
directed walks. 



Figure 2. (a) A directed walk, (b) A partially directed walk, (c) A spiral 
walk, (d) A prudent walk. 

The two simplest classes of SAW on the square lattice probably consist of directed and partially 
directed walks: a walk is directed if it involves at most two types of steps (for instance North and 
East), and partially directed if it involves at most three types of steps (Fig. [UJa-b)). Partially 
directed walks play a prominent role in the definition of our weakly directed walks. Among other 
solved classes, let us cite spiral SAW [371 HZ] an d prudent walks [U [THl [§]. We refer again to 
Fig. [5] for illustrations. Each time such a new class is defined and solved, one compares its 
properties to ((T|): have we reached with this class a large growth constant? Is the end-to-end 
distance of the walks sub-linear? 

At the moment, the largest growth constant (about 2.48) is obtained with prudent SAW. 
However, this is beaten by certain classes whose description involves a (small) integer k. like SAW 
confined to a strip of height k [TJ [35] , or SAW consisting of irreducible bridges of length at most 
k [201 12"5] . The structure of these walks is rather poor, which makes them rather unattractive 
from a combinatorial viewpoint. In the former case, they are described by a transfer matrix 
(the size of which increases exponentially with the height of the strip); in the latter case, the 
structure is even simpler, since these walks are just arbitrary sequences of irreducible bridges of 
small length. In both cases, the generating function is rational. The growth constant increases 
with k, providing better and better lower bounds on the growth constant of general SAW. The 
ability of solving these models for larger values of k mostly relies on progress in computer power. 
Regarding asymptotic properties, almost all solved classes of SAW exhibit a linear end-to-end 
distance, with the exception of spiral walks, which are designed so as to wind around their origin. 
But there are very few such walks [TTJ, as their growth constant is 1. 

With the weakly directed walks of this paper, we reach a growth constant of about 2.54. These 
walks are defined in the next section. Their generating function is given in Section [SJ after 
some preliminary results on partially directed bridges (Sections [3J and H]). This series turns out 
to be much more complicated that the generating functions of directed and partially directed 
walks, which are rational: we prove that it has a natural boundary in the complex plane, and 
in particular is not D-finite (that is, it does not satisfy any linear differential equation with 
polynomial coefficients). However, we are able to derive from this series certain asymptotic 
properties of weakly directed walks, like their growth constant and average end-to-end distance 
(which we find, unfortunately, to grow linearly with the length). Finally, we perform in Section[5] 
a similar study for a diagonal variant of weakly directed walks. Our intuition told us that this 
variant would give a larger growth constant, but we shall see that this is wrong. Section [7J 
discusses a few more points, including random generation. 

An extended abstract of this paper appeared in the proceedings of the 2010 FPSAC confer- 
ence [2]. 
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2. Weakly directed walks: definition 

Let us denote by N, E, S and W the four square lattice steps. All walks in this paper are 
self-avoiding, so that this precision will often be omitted. For any subset S of {N, E,S,W}, we 
say that a (self-avoiding) walk is an S-walk if all its steps lie in S. For instance, the first walk 
of Fig. [His a NE-walk, but also a NEW- walk. The second is a NEW- walk. We say that a SAW 
is directed if it involves at most two types of steps, and partially directed if it involves at most 
three types of steps. 

The definition of weakly directed walks stems from the following simple observations: 

(i) between two visits to any given horizontal line, a NE-walk only takes E steps, 

(ii) between two visits to any given horizontal line, a NEW- walk only takes E and W steps. 
Conversely, a walk satisfies (i) if and only if it is either a NE-walk or, symmetrically, a SE-walk. 
Similarly, a walk satisfies (ii) if and only if it is either a NEW- walk or, symmetrically, a SEW- 
walk. Conditions (i) and (ii) thus respectively characterize (up to symmetry) NE-walks and 
NEW- walks. 

Definition 1. A walk is weakly directed if, between two visits to any given horizontal line, the 
walk is partially directed (that is, avoids at least one of the steps N, E, S, \N). 

Examples are shown in Fig. [3] 



I 



J 



FIGURE 3. Two weakly directed walks. The second one is a bridge, formed of 5 
irreducible bridges. Observe that these irreducible bridges are partially directed. 

We will primarily focus on the enumeration of weakly directed bridges. As we shall see, this 
does not affect the growth constant. A self-avoiding walk starting at vo and ending at v n is a 
bridge if all its vertices v ^ v n satisfy h(vo) < h(v) < h(v n ), where h(v), the height of v, is 
its ordinate. Concatenating two bridges always gives a bridge. Conversely, every bridge can 
be uniquely factored into a sequence of irreducible bridges (nonempty bridges that cannot be 
written as the product of two nonempty bridges). This factorization is obtained by cutting the 
walk above each horizontal line of height n+ 1/2 (with n e Z) that the walk intersects only once 
(Fig. [21 right). It is known that the growth constant of bridges is the same as that of general 
self-avoiding walks [2S3 ■ The fact that bridges can be freely concatenated makes them useful 
objects in the study of self-avoiding walks [18l l20l l23l l24l [25] . 

The following result shows that the enumeration of weakly directed bridges boils down to the 
enumeration of (irreducible) partially directed bridges. It will be extended to general walks in 
Section [5) 

Proposition 2. A bridge is weakly directed if and only if each of its irreducible bridges is 
partially directed (that is, avoids at least one of the steps N, E, S, \N). In fact, this means that 
each of its irreducible bridges is a NES- or HSV\/-walk. 
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Proof. The second condition (being NES or NSW) looks more restrictive than the first one (being 
partially directed), but it is easy to see that they are actually equivalent: no non-empty ESW- 
walk is a bridge, and the only irreducible bridges among NEW- walks consist of a sequence of 
horizontal steps, followed by a N step: thus they are NES- or NSW- walks. 

So let us now consider a bridge whose irreducible bridges are partially directed. The portion of 
the walk lying between two visits to a given horizontal line is entirely contained in one irreducible 
bridge, and consequently, is partially directed. 

Conversely, consider a weakly directed bridge and one of its irreducible bridges w. Of course, 
w is also weakly directed. Let vo, ■ ■ ■ , v n be the vertices of w, and let Sj be the step that goes 
from Vi-i to Vi. We want to prove that w is a NES- or NSW-walk. Assume that, on the contrary, 
w contains a W step and an E step. By symmetry, we may assume that the first W occurs 
before the first E. Let Sk+i be the first E step, and let Sj be the last W step before Sk+i- Then 
Sj+i, . . . , Sfe is a sequence of N or S steps. Let h be the height of Sk+i- 

• Assume that Sj+i, • • • , Sfe are N steps (first walk in Fig. [J]). Let h! be the maximal height 
reached before Vj, say at Uj, with i < j. Then h' < h (otherwise, between the first 
visit to height h and Vk+i, the walk would not be partially directed). Given that w is 
irreducible, it must visit height hf again after Vk+i, say at V£. But then the walk joining 
Vi to ve is not partially directed, a contradiction. 

• Assume that Sj+i, • • • , Sfe are S steps (second walk in Fig. HJ. Let Vi, with i < k, be the 
last visit at height h before Vk- Then the portion of the walk joining Vi to Vk+i is not 
partially directed, a contradiction. 

Consequently, the irreducible bridge w is a NES- or NSW-walk. 




FIGURE 4. Illustrations for the proof of Proposition [5] 



We discuss in Section [5] a variant of weakly directed walks, where we constrain the walk to be 
partially directed between two visits to the same diagonal line (Fig. [5]). The notion of bridges 
is adapted accordingly, by defining the height of a vertex as the sum of its coordinates. We 
will refer to this model as the diagonal model, and to the original one as the horizontal model. 
There is, however, no simple counterpart of Proposition [5] a (diagonal) bridge whose irreducible 
bridges are partially directed is always weakly directed, but the converse is not true, as can be 
seen in Fig. [5] Thus bridges with partially directed irreducible bridges form a proper subclass of 
weakly directed bridges. We will enumerate this subclass, and study its asymptotic properties. 

3. Partially directed bridges: a step-by-step approach 

Let us equip the square lattice 1? with its standard coordinate system. With each model 
(horizontal or diagonal) is associated a notion of height: the height of a vertex v, denoted by 
h(v), is its ordinate in the horizontal model, and the sum of its coordinates in the diagonal 
model. Recall that a walk, starting at vq and ending at ti n , is a bridge if all its vertices v =/= v n 
satisfy h(vo) < h(v) < h(v n ). If the weaker inequality h(vo) < h(v) < h(v n ) holds for all v, we 
say the walk is a pseudo-bridge. Note that nonempty bridges are obtained by adding a step of 
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Figure 5 . Two weakly directed walks in the diagonal model. The second one is 
a bridge, factored into 6 irreducible bridges. Observe that the third irreducible 
bridge is not partially directed. 

height 1 to a pseudo-bridge (a N step in the horizontal model, a N or E step in the diagonal 
model). It is thus equivalent to count bridges or pseudo-bridges. 

By Proposition [21 the enumeration of weakly directed bridges in the horizontal model boils 
down to the enumeration of (irreducible) partially directed bridges. In this section and the 
following one, we address the enumeration of these building blocks, first in a rather systematic 
way based on a step-by-step construction, then in a more combinatorial way based on heaps of 
cycles. A third approach is briefly discussed in the final section. We also count partially directed 
bridges in the diagonal model, which will be useful in Section [51 




Figure 6. A N ES-pseudo-bridge in the horizontal model, (b) An ESW-pseudo- 
bridge in the diagonal model, (c) A N ES-pseudo-bridge in the diagonal model. 

As partially directed walks are defined by the avoidance of (at least) one step, there are four 
kinds of these. Hence, in principle, we should count, for each model (horizontal and diagonal), 
four families of partially directed bridges. However, in the horizontal model, there exists no 
non-empty ESW-bridge, and every NEW- walk is a pseudo-bridge. The latter class of walks is 
very easy to count, and has a rational generating function (Lemma If 6p . Moreover, a symmetry 
transforms NES-bridges into NSW-bridges, so that there is really one class of bridges that we 
need to count. In the diagonal model, we need to count ESW-bridges (which are equivalent to 
NSW-bridges by a diagonal symmetry) and NES-bridges (which are equivalent to NEW-bridges). 
Finally, to avoid certain ambiguities, we need to count ES-bridges, but this has already been 
done in [5]. 
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From now on, the starting point of our walks is always at height 0. The height of a walk is 
then defined to be the maximal height reached by its vertices. 

3.1. NES-BRIDGES IN THE HORIZONTAL MODEL 

Proposition 3. Let k > 0. In the horizontal model, the length generating function of NES- 
pseudo-bridges of height k is 



B^ k \t) 



t k 



G k (ty 

where Gk(t) is the sequence of polynomials defined by 

G_i = l, G = l-t, andfork>0, G k+l = (1 - t + t 2 + t 3 )G k - i 2 G fc _i. 
Equivalently, 

E v k t k _ ^ k r _ l-t- t 2 v 



or 



whe 



B (k) n\ z-^i - l-(l-t + t 2 + t 3 )v + t 2 v 2 

k>0 v ; fc>0 v ' 



((1 -t)U -t)U k - ({l-t)U-t)U k ' J 



T(t;u) = T(u) = J2m)i 



lT _ 1 - t + t 2 + t 3 - - t 4 ) (1 - 2t - t 2 ) 
~ It 

is a root of tu 2 — (l — t + t 2 + t 3 ) u + t = and U := 1/U is the other root of this polynomial. 

Proof. Fix k > 0. Let T be the set of NES- walks that end with an E step, and in which each 
vertex v satisfies < h(v) < k. Let 77 be the subset of T consisting of walks that end at height 
i. Let Ti(t) = Ti be the length generating function of %, and define the bivariate generating 
function 

k 

:>•■ 

i=0 

This series counts walks of T by their length and the height of their endpoint. Note that we 
often omit the variable t in our notation. The walks of Tk are obtained by adding an E step 
at the end of a pseudo-bridge of height k, and hence E>( k \t) = Th(t)/t. Alternatively, pseudo- 
bridges of height k containing at least one E step are obtained by adding a sequence of N steps 
of appropriate length to a walk of T, and this gives 

k 

B^(t)=t k + J2 T i( t ) tk ~ i =i fe (l + T(l/i)). (3) 

i=0 

(The term t k accounts for the walk formed of k consecutive N steps.) 

Lemma 4. The series T(t;u), denoted T(u) for short, satisfies the following equation: 

ut 2 t \ ml . 1 - (tu) k+1 (tu) k+1 , . , t 2 u 

1 T(u) = t — t T(l t) - T(t), 

1-tu 1 - tu J y ' 1-tu 1 - tu v ' ' 1 - tu w 

with u = 1/u. 

Proof. We partition the set T into three disjoint subsets, illustrated in Fig. [7] 

• The first subset consists of walks with a single E step. These walks read N • • • NE, with 
at most k occurrences of N, and their generating function is 

« l-{tu) k+1 

ty (tu) 1 = t — ^ — . 

i=0 
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• The second subset consists of walks in which the last E step is strictly higher than the 
previous one. Denoting by i the height of the next-to-last E step, the generating function 
of this subset reads 



k—i 



E,tu-(tu) k ut 2 , , (tu) k+1 m . 

i=0 \ j=l J i=0 



The third subset consists of walks in which the last E step is weakly lower than the 
previous one. Denoting by i the height of the next-to-last E step, the generating function 
of this subset reads 



1 — tu J 1 — tu 1 — tu 

i=0 \ j=0 / i=0 

Adding the three contributions gives the series T(u) and establishes the lemma. 





Figure 7. Recursive construction of bounded NES-walk in the horizontal model. 



The equation of Lemma 2] is easily solved using the kernel method (see e.g. [55]). The 
kernel of the equation is the coefficient of T(u), namely 

ut 2 t 
1 r- 

1 — tu 1 — tu 

It vanishes when u = U and u = U := 1/U, where U is defined in the lemma. Since T(u) is 
a polynomial in u, the series T{U) and T(C7) are well-defined. Replacing u by [/ or U in the 
functional equation cancels the left-hand side, and hence the right-hand side. One thus obtains 
two linear equations between T(t) and T(l/t): 



k+l 







t 



= i 



i - (tu) 

1 -tf/ 
1 - (tLQ fc+1 
1 -tU 



t 



(tU) 



k+l 



1-tU 



T(l/t) 



t 2 U 
1-tU 



(tU) k+1 t 2 U 
t^^T(l/t)--^-T(t). 



1-tU 



1-tU 



Solving this system gives in particular the value of T(l/t), and thus of B^(t) (thanks to ©). 
This provides the second expression of B^ k '(t) given in Proposition [31 The other results eas- 
ily follow, using standard connections between linear recurrence relations, their solutions, and 
rational generating functions [301 Thm. 4.1.1]. B 

3.2. ESW-BRIDGES IN THE DIAGONAL MODEL 

Proposition 5. Let k > 0. In the diagonal model, the length generating function of ESW- 
pseudo-bridges of height k is 



B[ k \t) 



t 



c k (ty 



where Gk(t) is the sequence of polynomials defined by 

G = l, G 1 = l-t 2 and for k > 1, G k+1 = (1 + t 2 )G k - t 2 {2 - t 2 )G fe _! . 
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Equivalently, 



\ - v k t k \ - , 1 - 2t 2 v 



or 



whe 



B?\t) = ^' (2 -^ - k 

(U - 2t)U k - ((2 - t 2 )U - 2t) ((2 - t 2 )U) 



u= l + t 2 -^(l-t 2 )(l-5t 2 ) 
2t 

is a root of tu 2 - (l + i 2 ) u + t (2 - t 2 ) = and U := 1/U. 

Proof. The proof is very close to the proof of Proposition [3l but the role that was played by E 
steps is now played by S steps. In particular, T is now the set of ESW-walks that end with a S 
step, and in which each vertex v satisfies < h(v) < k. The sets % and the series Tj(i) = Ti 
and T(t; u) = T(u) are then defined in terms of T as before. Note that Tk is in fact 0. 

Pseudo-bridges of height k containing at least one S step are obtained by adding a sequence 
of E steps of appropriate length to a walk of T, and © still holds (with replaced by B^). 

Lemma 6. The series T(t;u), denoted T(u) for short, satisfies the following equation: 

t 2 tu \ m/ , 9 1 - (tu) k , (tu) k , , . tu rr , l . 

1 T(u) = t 2 - t 2 - — —Til t) Tit), 

1-tu 1 - tu) ' 1-tu 1 - tu w ; 1 - tu w ' 

with u — 1/u. 

Proof. We partition the set T into three disjoint subsets. 

• The first subset consists of walks with a single S step. Their generating function is 

tuj-{t u y = t 2l - {tu)k . 

^ y ' 1-tu 

j=l 

• The second subset consists of walks in which the last S step is weakly higher than the 
previous one. Their generating function reads 



k—i 



tu £ | ItftV £(tu)J I = j^T(u) t 2 ^-T(l/t). 
i=o \ i=i / 

• The third subset consists of walks in which the last S step is strictly lower than the 
previous one. Their generating function reads 

i=0 \ j=0 J 

Adding the three contributions establishes the lemma. g 

Again, we solve the equation of Lemma [6] using the kernel method, and conclude the proof of 
Proposition using © (with B^ replaced by B[ k) ). B 

3.3. NES-BRIDGES IN THE DIAGONAL MODEL 

The generating function of NES-pseudo bridges is closely related to that of ESW-pseudo- 
bridges. This will be explained combinatorially in Section after a detour via partially directed 
excursions. 
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Proposition 7. Let k > 0. In the diagonal model, the length generating function of NES-pseudo- 

bridges of height k is 

{k) t^{2-ef 

* 2 (t) - G k (t) ' 
where Gk(t) is the sequence of polynomials defined in Proposition® In other words, 

Bi k \t) = (2-t 2 ) k B[ k \t), 

so that other expressions of B^ (t) can be derived from Proposition® 

Proof. The sets T, %, and the corresponding series T,(i) = Ti and T(t; u) = T{u) are defined as 
in the proof of Proposition [3] — only, the notion of height has changed. Note that To is 0. 

Pseudo-bridges of height k containing at least one E step are again obtained by adding a 
sequence of N steps of appropriate length to a walk of T, and still holds (with B^ replaced 
by ). 

Lemma 8. The series T(t\ u), denoted T{u) for short, satisfies the following equation: 

1 Ut ^ \t( )-t l ~( tu ^ k ( tu ^ k+1 T(l/t) ^ T(t) 

1 — tu 1 — tuJ 1 — tu 1 — tu 1 — tu ' 

with u = 1/u. 

Proof. We partition the set T into three disjoint subsets, defined as in the proof of Lemma |U 
• The first subset consists of walks with a single E step. Their generating function is 

k-i 



tu^2(tuy ^tu 1 



l-tu 

The second subset consists of walks in which the last E step is strictly higher than the 
previous one. Their generating function reads 

tu ^ ( my £ W ] = T^-T(u) - ^T(l/t). 



i=0 \ j=o 



• The third subset consists of walks in which the last E step is weakly lower than the 
previous one. Their generating function reads 

tu ]T U(t)u> = ^=T{u) - j^nt). 

1=0 \ j=i ) 

Adding the three contributions establishes the lemma. I 

The kernel occurring in Lemma [5] is obtained by replacing u by u in the kernel of Lemma [BJ 
Consequently, it vanishes when u = 1/U or u = U/(2 — i 2 ), where U is the series defined in 
Proposition [5] We use the kernel method to solve the equation of Lemma 13 and conclude the 
proof of Proposition [7] using © (with B^ replaced by B^). ■ 

3.4. ES-BRIDGES IN THE DIAGONAL MODEL 

We state our last result on partially directed bridges without proof, for two reasons. Firstly, 
the step-by-step approach used in the previous subsections should have become routine by now, 
and is especially simple to implement here. Secondly, this result already appears in [BJ Prop. 3.1] 
(where a bridge preceded by an E step is called a culminating path) . 
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Proposition 9. Let k > 0. In the diagonal model, the length generating function of ES-pseudo- 
bridges of height k is 

, > +k 

B " (<) ~ W 
where F k (t) is the sequence of polynomials defined by 

Equivalently, 



F-i = 1, F = 1, and for k > 0, F k+1 = F k - t 2 F k ^. 



\ - v k t k \ - , 1 -z;t 2 

^rOw"^" fc ~ 1 -v + v 2 t 2 ' 

k>0 n \ L ) k>0 



or 



where 



B ( k) (t) 



U 2 -U 2 

JJk+2 _ jjk+2 ' 



1 - x/T^4^ 

zs a root of tu 2 — u + 1 = and [7 := 1/fJ is t/ie other root of this polynomial. 

4. Partially directed bridges via heaps of cycles 

In this section, we give alternative (and more combinatorial) proofs of the results of Section [3] 
In particular, these proofs explain why the numerators of the rational series that count partially 
directed bridges of height k are so simple (t k or t k (2 — t 2 ) k , depending on the model). 

As a preliminary observation, let us note that ES-pseudo-bridges of height k in the diagonal 
model can be seen as arbitrary paths on the segment {0, 1, . . . , fc}, with steps ±1, going from 
to fc. Therefore, a natural way to count them is to use a classical result that expresses 
the generating function of paths with prescribed endpoints in a directed graph. This result is 
recalled in Proposition [TU] below. It gives a straightforward proof of Proposition [5) However, 
the other three classes of bridges that we have counted do not fall immediately in the scope 
of this general result, because of the self-avoidance condition (which holds automatically for 
ES- walks). For instance, in the horizontal model, a N ES-pseudo-bridge of height fc is not an 
arbitrary path with steps 0, ±1 going from to fc on the segment {0, 1, . . . , fc}. We show here 
how to recover the results of Section [3] by factoring bridges into more general steps, and then 
applying Proposition HU1 

Let r = (V, E) be a (finite) directed graph. To each arc of this graph, we associate a weight 
taken in some commutative ring (typically, a ring of formal power series) . A cycle of Y is a path 
ending at its starting point, taken up to a cyclic permutation. A path is self-avoiding if it does 
not visit the same vertex twice. A (non-empty) self-avoiding cycle is called an elementary cycle. 
Two paths are disjoint if their vertex sets are disjoint. The weight w(tt) of a path (or cycle) tt 
is the product of the weights of its arcs. A configuration of cycles 7 = {71, . . . ,7 r } is a set of 
pairwise disjoint elementary cycles. The signed weight of 7 is 

r 

«K7) :=(-l) r rh(7)- 

i=l 

For two vertices i and j, denote by Wij the generating function of paths going from from i to j: 

We assume that this sum is well-defined, which is always the case when Wij is a length generating 
function. 
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Proposition 10. The generating function of paths going from i to j in the weighted digraph T 
is 

where G = 10(7) is the signed generating function of configuration of cycles, and 

N i,j = ^2™(v)ti(l), 

■n,i 

where n is a self-avoiding path going from i to j and 7 is a configuration of cycles disjoint from 
V- 

This classical result can be proved as follows: one first identifies Wi j as the coefficient 
of the matrix (1 — A) -1 , where A is the weighted adjacency matrix of T. Thanks to standard 
linear algebra, this coefficient can be expressed in terms of the determinant of (1 — A) and one of 
its cofactors |30l Thms. 4.7.1 and 4.7.2]. A simple expansion of these as sums over permutations 
shows that the determinant is G, and the cofactor Ni j. Proposition [TU] can also be proved 
without any reference to linear algebra, using the theory of partially commutative monoids, or, 
more geometrically, heaps of pieces |15l I31| . In this context, configurations of cycles are called 
trivial heaps of cycles. This is the only justification of the title of this section, where no non-trivial 
heap will actually be seen. 

As a straightforward application, let us sketch a second proof of Proposition The vertices 
of r are 0,1, ... ,k, with an arc from i to j if \i — j\ = 1. We apply the above proposition to 
count paths going from to k. All arc weights are t. The elementary cycles have length 2, and 
by induction on k, it is easy to see that G, the signed generating function of configurations of 
cycles, is the polynomial F^. The only self-avoiding path r\ going from to k consists of k 'up' 
steps and visits all vertices, so that No,k = t k - Proposition [5] follows. 

4.1. Bridges with large down steps 

The proof of Proposition [9] that we have just sketched can be extended to paths with arbitrary 
large down steps. This will be used below to count partially directed bridges. 

Let Tfe be the graph with vertices {0, . . . , k} and with the following weighted arcs: 

• ascending arcs i — > i + 1 of height 1, with weight A, for i = 0, . . . , k — 1; 

• descending arcs i — > i — h of height h, with weight D^, for i = h, . . . , k and h > 0. 

For k > 0, denote by C^ k ' the generating function of paths from to k in the graph T^. These 
paths may be seen as pseudo-bridges of height k with general down steps. 

Lemma 11. The generating function of pseudo-bridges of height k is 

Ak 

where the generating function of the denominators Hk is 

1 - D(v, 
l-v + vD{vA) ' 



k>0 

with D{v) the generating function of descending steps: 

D(v) = D h v h . 

h>0 

Proof. With the notation of Proposition [TO1 the series C' fc - ) reads Nq^/G. Since all ascending 
arcs have height 1, the only self-avoiding path from to k consists of k ascending arcs, and has 
weight A k . As it visits every vertex of the graph, the only configuration of cycles disjoint from it 
is the empty configuration. Therefore, the numerator Nq^ is simply A k . The elementary cycles 
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consist of a descending step of height, say, h, followed by h ascending steps. The weight of this 
cycle is DhA h . 

To underline the dependence of our graph in k, denote by the denominator G of Proposi- 
tion [TU] Consider a configuration of cycles of IV either the vertex k is free, or it is occupied by 
a cycle; this gives the following recurrence relation, valid for k > 0: 

k 

Hk = Hk-i — ^ DhA h Hk-h-i, 

h=Q 

with the initial condition = 1. This is equivalent to ([4]). B 
4.2. Partially directed self-avoiding walks as arbitrary paths 

As discussed above, it is not straightforward to apply Proposition [TUJ (or Lemma [TTjl to the 
enumeration of partially directed bridges, because of the self-avoidance condition. To circumvent 
this difficulty, we will first prove that partially directed self-avoiding walks are arbitrary paths 
on a line with large down steps. 

It will be convenient to regard lattice walks as words on the alphabet {N, E, S, W}, and sets of 
walks as languages. We thus use some standard notation from the theory of formal languages [19] . 
The length of a word u (the number of letters) is denoted by and the number of occurrences 
of the letter a in u is \u\ a . For two languages C and CJ , 

• C + £' denotes the union of C and £'; 

• ££' denotes the language formed of all concatenations of a word of C with a word of 

• C* denotes the language formed of all sequences of words of C; 

• C + denotes the language formed of all nonempty sequences of words of C. 

Finally, for any letter a, we denote by a the elementary language {a}. A regular expression 
is any expression obtained from elementary languages using the sum, product, star and plus 
operators. It is unambiguous if every word of the corresponding language has a unique factor- 
ization compatible with the expression. To take a simple example, the expressions (N + E)* and 
(N + W)* are unambiguous expressions describing NE- and NW- walks respectively. However, the 
expression (N + E)* + (N + W)* is ambiguous, as every N-walk is matched twice. Unambiguous 
regular expressions translate directly into enumerative results. 

Let us say that a NES-walk is proper if it neither begins nor ends with a S step. All NES- 
pseudo-bridges are proper, whether in the horizontal or diagonal model. The following lemma 
explains how to see proper NES-walks as sequences of generalized steps. 

Lemma 12. Every proper NES-walk has a unique factorization into N steps and nonempty 
proper ES-walks with no consecutive E steps. In other words, the language of proper NES-walks 
admits the following unambiguous regular expression: 

(N + E(S+E)*)*. 

Proof. The factorization of proper NES-walks is exemplified in Fig. [HJ Every N step is a factor, 
as well as every maximal ES-walk with no consecutive E steps. B 



-• — •- 



Figure 8. The factorization of a proper NES-walk. 
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A similar result holds for ESW- walks (which we need to study in the diagonal model), which 
are obtained by applying a quarter turn to NES-walks. Let us say that an ESW- walk is proper 
if it neither begins nor ends with a W step. After a rotation, Lemma [T2] gives for the language 
of proper ESW-walks the following unambiguous description: 

(E + S(W+S)*)*. (5) 
4.3. Partially directed bridges 

We can now give new proofs of the results of Section [3J based on Lemma ITT1 

Second proof of Proposition® Thanks to Lemma [121 self- avoiding N ES-pseudo-bridges of height 
k can be seen as arbitrary pseudo-bridges of height k (in the sense of Lemma ITTj) where N is the 
only ascending step (of height 1 and weight i), and all words of E(S + E)* are descending steps. 
Moreover, the weight of a descending step u is t> u > and its height is \u\s- Thus, with the notation 
of Lemma 111! A = t and the generating function D(v) of descending steps is derived from the 
regular expression E(S + E)*: 

D(v) = ^5-. 

1 l-tv 

Proposition (3J in the form ([2]). now follows from Lemma [TT] g 

Second proof of Proposition^ Thanks to ([5]), self- avoiding ESW-pseudo-bridges of height k can 
be seen as arbitrary pseudo-bridges of height k (in the sense of Lemma fTTj) where E is the only 
ascending step (of weight t), and all words of S(W + S)* are descending steps. Moreover, the 
weight of a descending step u is t> u > and its height is |«|. Thus, with the notation of Lemma [TT1 
A = t and the generating function D(v) of descending steps is 

D(v) = — * 

1 l-tv 

Proposition [5] now follows from Lemma [TT] g 

Second proof of Proposition [7J Again, the description of N ES- walks given by Lemma [T2] allows 
us to regard these self-avoiding walks as arbitrary paths with generalized steps. In the diagonal 
framework, the ascending steps u are N and all words of E(SE)*. They all have weight t> u >. All 
other words of E(S + E)* are descending. Moreover, the weight of a descending step u is t'"l and 
its height is |w|s — |u|e- Thus, with the notation of Lemma [TT1 

t t(2-t 2 ) , . tv- 1 tv- 1 

A = (+ r? = W and d{v) = y^-^- 

l-tv 

However, one must pay attention to the following detail: in a N ES-pseudo-bridge of height 
k, only the last generalized step can end at height k, because all descending steps begin with 
E. Similarly, all descending steps end with E, which implies that the only generalized step that 
starts at height is the first one (and moreover it is an ascending step). Thus a N ES-pseudo- 
bridge of height k > 2 is really a pseudo-bridge (in the sense of Lemma (TTJ of height k — 2, 
preceded and followed by an ascending step. Thus for k > 2, 

D 2 



where the generating function of the denominators is given in Lemma [TT1 Given that Bq 



Hk-2 

(2) 



f 2) 

and B\ J — A, we have 



H2~t 2 ) k v k t(2~t 2 )v T-^t k (2-t 2 ) k v k l-2t 2 v 

~ 1 + X -4-^ + E A* Hk - 2 - 



B f) {t) A ^ A k l-(l + t 2 )v + t 2 (2-t 2 )v 2 

This is equivalent to Proposition [7] 
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5. Weakly directed walks: the horizontal model 

We now return to the weakly directed walks defined in Section [2] We determine their gen- 
erating function, study their asymptotic number and average end-to-end distance, and finally 
prove that the generating function we have obtained has infinitely many singularities, and hence, 
cannot be D-finite. 

5.1. Generating functions 

By combining Propositions [5] and [3J it is now easy to count weakly directed bridges. 
Proposition 13. In the horizontal model, the generating function of weakly directed bridges is: 

wit) = TT7 ^ nr 



1+tB 

where B := ^2 k>0 B^ k '(t) is the generating function of NES-pseudo-bridges, given by Proposi- 
tion^ 

Proof. Let 2e be the set of irreducible NES-bridges, and let /e^) be the associated length 
generating function. We will most of the time omit the variable t in our series, writing for 
instance Je instead of /e(0- Given that a non-empty NES-bridge is obtained by adding a N step 
at the end of a NES-pseudo- bridge, and is a (non-empty) sequence of irreducible NES-bridges, 
we have: 

*B =-*=-. 

Define similarly the set Zyv, and the associated series 7w- By symmetry, Jvv = Ie- Moreover, 

i E ni w = N. 

Hence the generating function of irreducible bridges that are either NES or NSW is 

I:=I E + I W -t = ^L-t. 

By Proposition the generating function of weakly directed bridges is W = j— r . The result 
follows. I 

We will now determine the generating function of (general) weakly directed walks. As we 
did for bridges, we factor them into "irreducible" factors, but the first and last factors are not 
necessarily bridges, so that we need to extend the notion of irreducibility to more general walks. 
Let us say that a walk vq ■ ■ ■ v n is positive if all its vertices v satisfy h(v) > h(vo), and that it is 
copositive if all vertices v ^ v n satisfy h(v) < h{v n ). Thus a bridge is a positive and copositive 
walk. 

Definition 14. Let r denote the reflection through the x-axis. A non-empty walk w is N- 
reducible if it is of the form qp, where q is a nonempty copositive walk and p is a nonempty 
positive walk. It is S-reducible if r(w) is N-reducible. Finally, it is irreducible if it is neither 
N-reducible nor S-reducible. 

We can rephrase this definition as follows. If a horizontal line at height h + 1/2, with /igZ, 
meets w at exactly one point, we say that the step of w containing this point is a separating 
step. Of course, this step is either N or S. Then a non-empty walk is irreducible if it does not 
contain any non-final separating step. It is then clear that the above definition extends the 
notion of irreducible bridges defined in Section [5] a non-empty bridge is never S-reducible, and 
it is N-reducible if and only if it is the product of two non-empty bridges. Also, observe that the 
endpoint of a N-reducible walk is strictly higher than its origin: Thus a walk may not be both 
N-reducible and S-reducible. 

By cutting a walk after each separating step, one obtains a decomposition into a sequence of 
irreducible walks. This may be either a N-decomposition or a S-decomposition. The first factor 
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of a N-decomposition is copositive, while the last one is positive. The intermediate factors are 
bridges. 

We can now generalize Proposition [3J and characterize weakly directed walks in terms their 
irreducible factors. 

Proposition 15. A walk is weakly directed if and only if each of its irreducible factors is partially 
directed. Equivalently, each of these factors is a NES- or a HS\N-walk. 

Proof. The proof is very similar to that of Proposition^ First, the equivalence between the two 
conditions comes from the fact that any partially directed irreducible walk is NES or NSW. 

Now, if all irreducible factors of a walk are partially directed, then this walk is weakly directed: 
two points of the walk lying on the same horizontal line belong to the same irreducible factor. 

Conversely, let w be an irreducible factor of a weakly directed walk; then w is weakly directed. 
We prove that it is either a NES- or a NSW-walk. Assume that this is not the case, i.e. w contains 
both a W and an E step. By symmetry, we may assume that it contains a W step before its 
first E step. By symmetry again, we may assume that, between the first E step and the last W 
step that precedes it, the walk consists of N steps. Then the first argument used in the proof of 
Proposition [2j depicted in the first part of Fig. |4l leads to a contradiction. g 

We now proceed to the enumeration of general weakly directed walks. 

Lemma 16. The generating functions T{t), P(t) and Q(t) of general, positive and copositive 
NES-walks are: 

l + t 



T(t) 



l-2t-t 2 ' 




1 -t* 

-1-t 



\-2t-t 2 
Q(t) = l+tP{t). 

Proof. Let us start with general NES-walks. The language T of these walks is given by the 
following unambiguous description: 

T=N*+S+ + TE(N*+S+), 

from which the expression of T(t) readily follows. 

Let us now count positive walks. Let P(t; u) be their generating function, with the variable u 
accounting for the height of the endpoint. We decompose positive walks by cutting them before 
the last E step; this is similar to what we did in the proof of Lemma HJ We thus obtain: 

pu. u)= 1 , t 2 uP(t;u) | t(P(t;u)-tuP(t;t)) 
' 1 — tu 1 — tu 1 — tu 

We rewrite this as follows: 

U P(t;lt) - 1 ^n*,t) 



1 — tu 1 — tu J 1 — tu 1 — tu 

We apply again the kernel method: we specialize u to the series U of Proposition [3J this cancels 
the coefficient of P(t;u), and we thus obtain the value of P(t;t). We then specialize the above 
equation to u = 1 to determine P(t; 1), which is the series denoted P(t) in the lemma. 

Finally, a non-empty copositive walk is obtained by reading a positive walk, seen as a word 
on {N, E, S, W}, from right to left, and adding a final N step. This gives the last equation of the 
lemma. B 

Proposition 17. The generating function of weakly directed walks is 

W{t) = 1 + (2T<(t) - 2b) + 2(2Q<(t) - t)W(t)(2P(t) - t), 
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where the series Ti(t), Pi{t) and Qi{t) count respectively general, ■positive, and copositive irre- 
ducible NES-walks, and are given by: 

Ti(t) = T(t) 1 - 2Q i (t)(l + tB(t))P, t (t), P t (t) = ^l B{ l y Qi(t) = YTWWY 
The series W, B, T, P and Q are those of Proposition ^^ and Lemma \W\ 

Proof. In order to determine the series Ti(t), Pi(t) and Qi(t), we decompose into irreducible 
factors the corresponding families of NES-walks. 

• A general NES-walk is either 

— empty, or 

— irreducible, or 

— N-reducible: in this case, it consists of an irreducible copositive NES-walk, followed 
by a sequence of irreducible N ES-bridges (forming a N ES-bridge) , and finally by an 
irreducible positive NES-walk; or 

— symmetrically S-reducible. 

Since 1 + tB(t) is the generating function of bridges, this gives 

T(t) = 1 + Tt(t) + 2Qi(t)(l + tB{t))Pi{t). 

• We now specialize the above decomposition to positive NES-walks. Observe that when 
such a walk is N-reducible, its first factor is a bridge. This allows us to merge the second 
and third cases above. Moreover, the fourth case never occurs. Thus a positive NES-walk 
is either 

— empty, or 

— a N ES-bridge followed by an irreducible positive NES-walk. 
This yields: 

P(t) =l+(l + tB(t))Pi(t). 

• We proceed similarly for copositive NES-walks. Such a walk is either 

— empty, or 

— an irreducible copositive NES-walk followed by a N ES-bridge. 
This gives: 

Q(t) = l + Qi(t)(l + tB(t)). 
We thus obtain the expressions of Ti, Pj and Qi announced in the proposition. 

Recall from Proposition [T5] that a walk is weakly directed if and only if its irreducible factors 
are NES- or NSW- walks. Thus, 

• a weakly directed walk is either 

— empty, or 

— an irreducible NES- or NSW- walk, or 

— N-reducible: it then factors into an irreducible copositive NES- or NSW- walk, a 
sequence of NES- or NSW- irreducible bridges (forming a weakly directed bridge), 
and an irreducible positive NES- or NSW- walk; or 

— symmetrically S-reducible. 

The contribution to W(t) of the first case is obviously 1. The only irreducible walks that are 
both NES and NSW are N and S. The generating function of irreducible NES- or NSW- walks 
is thus 2Ti(t) — 2t. Similarly, the generating function of irreducible positive (resp. copositive) 
NES- or NSW-walks is 2Pi(t) — t (resp. 2Qi(t) — t). In both cases, the term —t corresponds to 
the walk reduced to a N step, which is both NES and NSW. Adding the contributions of the four 
classes yields the announced expression of W(t). a 
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5.2. Asymptotic results 

Proposition 18. The generating function W of weakly directed bridges, given in Provosition lllil 
is meromorphic in the disk T> = {z : \z\ < y/2 — 1}. It has a unique dominant pole in this disk, 
p ~ 0.3929. This pole is simple. Consequently, the number w n of weakly directed bridges of 
length n satisfies 

w„ ~ Kp n , 

with fj, = l/p~ 2.5447. 

Let N n denote the number of irreducible factors in a random weakly directed bridge of length 
n. The mean and variance of N„ satisfy: 

E(JV n )~mn, Y(N n )~s 2 n, 

where 

m~ 0.318 and s 2 ~ 0.7, 

and the random variable N " ! ~^ n converges in law to a standard normal distribution. In particu- 
lar, the average end-to-end distance, being bounded from below by K(N n ), grows linearly with n. 
These results hold as well for general weakly directed walks, with other values of n, m and 5. 

Proof. Recall from the proof of Proposition IT51 that Wit) = 1/(1 — I it)), where I(t) counts 
partially directed irreducible bridges, which are certain NES- or NSW- walks. The generating 
function T(t) of NES-walks, given in Lemma [TBI has radius of convergence s/2— 1. Hence, I has 
radius of convergence at least v2 — 1, and W is meromorphic in the disk V. 

In this disk, we find a pole at each value of t for which I(t) = 1. As lit) has non- negative 
coefficients and is aperiodic, a pole of minimal modulus, if it exists, can only be real, positive 
and simple. Thus if there is a pole in T>, then W has a unique dominant pole p, which is simple, 
and the asymptotic behaviour of the numbers w n follows. 

In order to prove the existence of p, we use upper and lower bounds on the series lit). 
For any series F(t) = Y, m > f m t m , and n > 0, denote F< n (t) := £m=o /™* m and F >n(t) :== 
E m >„ fmt m - Then for < t < V2 - 1 and n > 0, we have 

I-(t)<Iit)<I + it), (6) 

where the series 

l~(t) := I< n (t) and J+(t) := /<„(*) + 2T >n (t) = /<„(*) + 2T(t) - 2T<„(t) 

can be evaluated exactly for a given value of n. The upper bound follows from the fact that I 
counts walks that are either NES- or NSW- walks. Using these bounds, we can prove the existence 
of p and locate it. More precisely, 

P~ < P < P + , (7) 

where 

/-(p+) = i+0o = i. 

Taking n = 300 gives 5 exact digits in p, = 1/p. 

Let us now study the number of irreducible bridges in a random weakly directed bridge. The 
series that counts these bridges by their length and the number of irreducible bridges is 

One easily checks that Wit, x) corresponds to a supercritical sequence, so that Prop. IX. 7 of |14| 
applies and establishes the existence of a gaussian limit law, after standardization. Regarding 
the estimates of m and s, we have 

1 and + 
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As I(t) has non-negative coefficients, we can combine the bounds ((SJ) on I(t) and ([7]) on p to 
obtain bounds on the values of m and s. 

Consider now the generating function W of general weakly directed walks, given in Proposi- 
tion [TTJ The series X^, Pi and Qi count certain partially directed walks, and thus have radius at 
least a/2-1. Moreover, 2Q t (t) - t> and 2Pi(t) - t > for t > 0. Hence W has, as W itself, 
a unique dominant pole in T>, which is p. 

The argument used to prove Proposition [T7] shows that the series that counts weakly directed 
walks by their length and the number of irreducible factors is 

W(t,x) = l + x(2T l (t) - 2t) + 2x 2 {2Qi(t) - t)W(t,x){2Pi{t) -t), 

where W(t,x) is given by ([5]). This yields the announced results on the number of irreducible 
factors. a 

5.3. Nature of the series 

Proposition 19. The generating function B = ^2 k>0 B^ (t) of NES-pseudo-bridges, given in 
Proposition^ converges aroundO and has a meromorphic continuation inC\£, where £ consists 
of the two real intervals [—a/2 — 1,-1] and [a/2 — 1,1], and of the curve 

f • „ 2 l-x 2 -2x 3 } 

£o=^ + i y:x>0, y = 1 + 2x )■ 

This curve, shown in Fig. [9j is a natural boundary of B . That is, every point of£o is a singularity 
ofB. 

The above statements hold as well for the generating function W of weakly directed bridges, 
given in Proposition [T51 In particular, neither B nor W is D-finite. 
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FIGURE 9. The curve £q and the zeroes of G*2o. 

Before proving this proposition, let us establish two lemmas, dealing respectively with the 
series U and the polynomials Gk occurring in the expression of B (Proposition [3]) . 

Lemma 20. For t € C \ {0}, the equation t(u + 1/u) — 1 — t + t 2 + t 3 has two roots, counted 
with multiplicity. The product of these roots is 1. Their modulus is 1 if and only if t belongs to 
the set £ defined in Provosition \19l 
Let 

l-t + t 2 +t 3 -^(l-t+)(l-2t-t2) 
[) ~ 2t 
be the root that is defined at t = 0. This series has radius of convergence a/2-1. It lias 
singularities at ±\/2— 1, ±1 and ±i, and admits an analytic continuation in 

C \ ([-V2 - 1, -1] U [a/2 - 1, 1] U [i, ioo) U [-i, -ioo)) . 

Proof. The first two statements are obvious. Now assume that the roots u and 1/u have modu- 
lus 1, that is, u = e lB for 6 £ R. This means that f(t) := 1 ^ t "^* +t = cosO is real, and belongs 
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to the interval [—1, 1]. Write t = x + iy, and express Im f(t) in terms of x and y. One finds that 
f(t) is real if and only if either y = (that is, t G K) or 

y 2 (l + 2x) = l~x 2 -2x 3 . (9) 

Since y 2 > 0, this is only possible if —1/2 < x < x c where x c ~ 0.65... satisfies 1 — x 2 — 2 x 3 = 0. 
Observe that the above curve includes £$. 

For real values of t, an elementary study of / shows that f(t) G [—1, 1] if and only if t G 
[-y/2 - 1,-1] U [y/2 - 1, 1] (see Fig. QUI left). If t = x + iy is non-real and © holds, then 
f(t) = ~ 1+ 1 4 j c ~+ 4x . Given that -1/2 < x < x c , this belongs to [—1,1] if and only if x is 
non- negative (see Fig. [TU1 middle). We have thus proved that |u| = 1 if and only if t G £. 

The properties of the series U follow from basic complex analysis. Of course, one may choose 
the position of the cuts differently, provided they include the 6 singularities. With the cuts along 
the coordinate plot of the modulus of U is shown on the right of Fig. [TUJ m 




FIGURE 10. The functions t ^ f(t) = 1 ~ t+ ^ +t \ x ^ ^{f^^ ' and a P lot 
of the modulus of U, showing the two cuts on the real axis. 



Lemma 21. Let £ be the subset of C defined in Proposition 1191 and Gk the polynomials of 
Proposition [31 

IfGk(z) = Gt(z) = with t ^ k, then z e £. 
If Gk(z) — and z is non-real, then z ^ £. 

The set of accumulation points of roots of the polynomials Gk is exactly £ . 

The latter point is illustrated in Fig. O 
Proof. Note first that, for z ^ 0, 



G k (z) 



((1 — z)u — z)u k — ((1 — z)/u — z) 



-k 



u — 1/u 

where u and 1/u are the two roots of z(u + 1/u) = 1 — z + z 2 + z 3 . 

Assume Gk{z) = Gi(z) = 0. Then z ^ (because Gk(0) = 1). The equations Gk(z) = 
Ge(z) = imply u 2k = u 2e , so that |u| = 1, that is, by Lemma |2"01 z G £. 

Assume Gk(z) = and z is non-real. Assume moreover that z G £ . Let u and 1/u be defined 
as above. By Lemma I2TJ1 |ti| = 1. Write u = e lS . Then Gk(z) = implies 

z _ sin((fc+l)fl) 
1 — 2 sm(k6) 
which contradicts the assumption that z is non-real. 
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Now let 2 be an accumulation point of roots of the Gfc's. There exists a sequence Zj that 
tends to z such that G^izi) = 0, with fej — > oo. We want to prove that z E £. If z is one of the 
6 singularities of U, then there is nothing to prove. Otherwise, U has an analytic description in 
a neighborhood of z. The equation Gk f {zi) = reads 

rj, «,2k, _ (1 - Zi)/U{zi) - g 

By continuity, Z7(zi) — ► as i — > oo. If z = 0, then U(z) = and the right-hand side diverges 
while the left-hand side tends to 0. This is impossible, and hence z / 0. This implies that 
the right-hand side tends to a finite, non-zero limit and, by continuity, forces |f(z)| = 1. By 
Lemma 1201 this means that z £ £. 

Conversely, let z £ £. By Lemma [201 the two roots of z(u+l/u) = l — z+z 2 +z 3 can be written 
e ±i8 _ gy density, we may assume that 9 = jit/d, for < j < I. This excludes in particular the 
6 singular points of U , for which u = ±1. This means that U has an analytic description in a 
neighborhood of z, so that for t close to z, U(t) = U(z)(l + s) with s = (t—z) jjnrj + 0((t — z) 2 ). 
Thanks to the equation satisfied by U(z), it is easy to see that U'(z) ^ if z ^ x c , where x c is 
defined in the proof of Lemma 1201 We assume from now on that z x c (again, by density, this 
is a harmless assumption). Let k be a multiple of L The equation Gk(t) = reads 

r T(A 2k _ (l-t)/U(t)-t 

W (l-t)U(t)-t> 

or, given that U(z) 2k = e 2t: > k7T ^ t = 1, 

(l + s)"= { }- z) (™- Z + 0(t-z). 
(1 — z)U (z) — z 

The right-hand side being finite and non-zero, one finds a root t of Gk in the neighborhood of z: 

U(z) , f{l-z)/U(z)-z\ . 

and this root gets closer and closer to z as k increases. Thus z is an accumulation point of roots 
of the Gk's. ■ 

Proof of Proposition 1191 

One has = ^t fc /G fe (t), with 

< fc it — 1 /u 



G fe (i) ((1 - t)u - t)u k - ((1 - t)/u - t)i 



(10) 



u and 1/u being the roots of t(u + 1/it) = 1 — t + t 2 + t 3 . Let us first prove that this series 
defines a meromorphic function in C\£. Assume t g 1 £ . By Lemma [zTTl t is not an accumulation 
point of roots of the polynomials Gk, and cancels at most one of these polynomials. Hence 
there exists a neighborhood of t in which at most one of the G^'s has a zero, which is t itself. 
Moreover, by Lemma [201 one of the roots u and 1/u has modulus larger than 1. By continuity, 
this holds in a (possibly smaller) neighborhood of t. Then (|10[) shows that the series t k /Gk(t) 
is meromorphic in the vicinity of t. Given that C \ £ is connected, we have proved that this 
series defines a meromorphic function in C \ £. The same holds for W(t), which is a rational 
function of t and B{t). 

Let us now prove that £q is a natural boundary of B. By Lemma [2TJ every non-real zero of 
Gk is in C \ £, and does not cancel any other polynomial Gi. Hence it is a pole of B. Now let 
z G £q with z $l M. By Lemma 1211 this point is an accumulation point of (non-real) zeroes of 
the polynomials Gk, and thus an accumulation point of poles of B. Thus it is a singularity of 
B, and the whole curve £q is a natural boundary of B. Given that B and W are related by a 
simple homography, this curve is also a boundary for W. ■ 
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6. The diagonal model 

We have defined weakly directed walks in the diagonal model by requiring that the portion 
of the walk joining two visits to the same diagonal is partially directed. This is analogous to 
the definition we had in the horizontal model. The definition of bridges is adapted accordingly, 
by defining the height of a vertex as the sum of its coordinates. However, there is no simple 
counterpart of Proposition [2j the irreducible bridges of a weakly directed bridge may not be 
partially directed (Fig. [5]). However, it is easy to see that bridges formed of partially directed 
irreducible bridges are always weakly directed. In this section, we enumerate these walks and 
study their asymptotic properties. 

6.1. Generating function 

Proposition 22. The generating function of bridges formed of partially directed irreducible 
bridges is 

Wa ^ = " 2tB[ WB 2 2tB Q ' 



1 + 2i — 



1 + £Bi 1 + 2tB 2 1 + tB a 



where the series Bi = ^2 k>0 B^ k \t) are given in Propositions [SJ [7J O 

Proof. Let Is be the set of irreducible ESW-bridges, and let Is be the associated length generating 
function. Given that a non-empty ESW-bridge is obtained by adding an E step at the end of a 
ESW-pseudo-bridge, and a non-empty sequence of irreducible ESW-bridges, there holds 

Define similarly the sets In, I e and Iw, and the associated series In, 7 e and iyv- Finally, let 2es 
(resp. Xnw) be the set of irreducible ES-bridges (resp. NW-bridges) , and let Ies (resp. Inw) be 
the associated series. Then 

2tB 2 = lE T and tB 



l-h l-hs 
(The factor 2 comes from the fact that a NES-bridge may end with a l\l or E step.) By symmetry, 
In = Ie, Iw = Is and Xes = X\iw- Moreover, 

i s nx N = e, i s nx E = Xes, Xs nx w = 0, x E nx N = n + e, x e nx w = n, x w nx N = x NW . 

By an elementary inclusion-exclusion argument, the generating function of partially directed 
irreducible bridges is 

2tB 1 AtB 2 2tB Q 

I := 2I S + 2I E - 2I ES -2t= — + 2t. 

b b 1 + tB x 1 + 2tB 2 1 + tB 

Hence the generating function of bridges formed of partially directed irreducible bridges is Wa = 
j^j. The result follows. I 

6.2. Asymptotic properties 

We obtain for the diagonal model asymptotic results that are similar to those obtained in the 
horizontal model, with a slightly smaller growth constant. We have to confess that this contra- 
dicts our original intuition: since in the horizontal model, two of the four classes of irreducible 
partially directed bridges (namely, ESW and NEW) are either trivial or degenerate, while in the 
diagonal model, all four classes are non-trivial, we thought we had a chance to reach a better 
growth constant in the diagonal model. This is unfortunately not the case. We nonetheless 
present this diagonal variant, because we believe it to be a natural attempt. We analyze below 
what makes the difference between the two growth constants, and this analysis shows that our 
hopes could just as well have come true. 



22 



A. BACHER AND M. BOUSQUET-MELOU 



Proposition 23. The generating function W/\ given by Proposition [55] is meromorphic in the 
disk V = {\z\ < y/2- 1}. It has a unique dominant pole in this disk, at p\ ~ 0.3940. This pole 
is simple. Consequently, the number of n-step bridges formed of partially directed irreducible 
bridges is asymptotically equivalent to K \x n , with fi = 1/pi ~ 2.5378. 

Let N n denote the number of irreducible bridges in a random n-step bridge formed of partially 
directed irreducible bridges. The mean and variance of N n satisfy: 

E(iV n )~mn, V(iV„)~s 2 n, 

where 

m~ 0.395 and s 2 = l±2.10~ 3 , 

and the random variable Nn B ^ n converges in law to a standard normal distribution. In particu- 
lar, the average end-to-end distance, being bounded from below by E(N n ), grows linearly with n. 

Proof. The arguments are the same as in the proof of Proposition [T5J This series reads Wa = 
1/(1 — 7), where I counts partially directed irreducible bridges. The only change is in the bounds 
we use on the series /: 

I-(t)<I(t)<I + (t), 

with 

r(t) := 7< n (t) and I+(t) := /<„(<) + 4T >n (t) = /<„(*) + 4T(t) - 4T< n (i), 
where T(t) is the generating function of NES-walks, given in Lemma [Tin g 

Remark. Hence the growth constant in the diagonal model is a bit smaller than in the horizontal 
model. This does not seem to be predictible. The series of Propositions IT31 and l2"2l respectively 
read 

w ®=T^m and ^ (<) = i—b) 

where I and I a count irreducible partially directed bridges, respectively in the horizontal and 
diagonal model. As t increases from to s/2— 1 (the radius of convergence of the series of partially 
directed walks), 7a(£) = 2£ + 0{t 2 ) first dominates I(t) = t + 0(t 2 ), but the graphs of these two 
functions cross before any of them reaches 1 (Fig. [TTjl . so that I(t) reaches 1 before Ja(0 does. 
The fact that the graphs cross is consistent with our belief that / has radius 1 ~ 0.41 while 
I A has a larger radius of convergence, namely l/y/E ~ 0.44. But /a(0 could just as well have 
reached 1 before the crossing point. 



FIGURE 11. The functions I(t) and I A (t) for t G (0.35, 0.4). The function I A (t) 
first dominates, but the graphs cross before the functions reach 1. 
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7. Final comments 

7.1. One more way to count partially directed bridges 

We have presented in Sections [3] and 2] two approaches to count partially directed bridges. 
Here, we discuss a third method, based on standard decompositions of lattice paths. This 
approach involves some guessing, whereas the two others don't. In the diagonal model, it allows 
us to understand more combinatorially why the generating functions of NES- and ESW-bridges 
only differ by a factor (2 — t 2 ) k . Moreover, this approach is needed for random generation. 

We first discuss the horizontal model, that is, the enumeration of NES-bridges given by Propo- 
sition [21 We say that a NES-walk is an excursion if it starts and ends at height 0, and all its 
vertices lie at a non-negative height. Let E^ = E^ k \t) be the length generating function of 
excursions of height at most k. As before, B^ k ' denotes the generating function of NES-pseudo- 
bridges of height k. 




FIGURE 12. Recursive decomposition of N ES-excursions and pseudo-bridges. 

Excursions and pseudo-bridges can be factored in a standard way by cutting them at their 
first (resp. last) visit at height 0. These factorizations are schematized in Fig. [T2] They give: 

• for excursions of bounded height, the recurrence relation 

= 1 + tE™ + t 2 - l) + 1 3 - l) , 

with the initial condition = 1, 

• for pseudo-bridges of height fc, the recurrence relation 

B (k) = (i + tE^tB^-V, 

with the initial condition B^ = 1/(1 - t). 
It is then straightforward to check by induction on k that 

where Gk is the sequence of polynomials defined in Proposition [3J Of course, these expressions 
have to be guessed — which is actually not very difficult using a computer algebra system. 

Let us now discuss the diagonal model. We will recover Proposition [5] (for ESW-bridges, 
or, equivalently, NSW-bridges) and Proposition [7] (for NES-bridges). Let E^ and E^ be the 
generating functions of NSW- and N ES-excursions, respectively, of height at most fc. These two 
series coincide: indeed, a NES-excursion is obtained by reading a NSW-excursion backwards, 
reversing each step, and this does not change the height of the excursion. 

Partially directed excursions and pseudo-bridges can be factored as in the horizontal 
illustrated in Fig. [121 This gives: 

• for NSW-excursions, the recurrence relation 

E[ k) = l + t 2 (E[ k - 1] - 1) + t 2 E[ k - 1] E{ k ) , 
with the initial condition E^ = 1, 
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FIGURE 13. Top: Recursive decomposition of NSW-excursions (or equivalently, 
of NES-excursions). Bottom: Recursive decomposition of NSW-pseudo-bridges 
(left) and NES-pseudo-bridges (right). 

• for NSW-pseudo-bridges, 

B (k) = D (k) tB (k-i) ^ (n) 

where counts NSW-excursions of height at most k not ending with a S step, and 

B? ] = 1, 

• for NES-pseudo-bridges, 

B { 2 k) = (l + E^tB^, (12) 

with the initial condition B^ = 1. 
As in the horizontal model, it is easy to check by induction that 

(2 -t 2 )%I-l and (2 ~' 2) *' fc , 

where Gk is the sequence of polynomials defined in Proposition [5j This gives Proposition [7] 
In order to prove Proposition [5j we will prove combinatorially that 1 + = (2 — t 2 )D^\ 

By comparing {TTJ and fT2]l. this will establish the link B { 2 h) = (2 - t 2 ) k B[ h) between the two 

types of pseudo-bridges. 

First, let u be a NSW-excursion of height at most k, ending with NW. Writing u = uNW, we 

see that v is an excursion that does not end with a S step; therefore, excursions ending with NW 

are counted by t 2 D± '. 

Let now u be an arbitrary NSW-excursion of height at most k. We distinguish two cases: 

• either u does not end with a S step; such excursions are counted by D^; 
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• or u reads vS; then v does not end with a N step. Let v! = iA/V: then u' is an excursion 
ending with W but not with NW. According to the above remark, such excursions are 
counted by D[ k) - 1 - t 2 D [k) . 

Putting this together, we find = (2 — t 2 )D^ — 1, which concludes the proof. 



7.2. Random generation of weakly directed bridges 

We now present an algorithm for the random generation of weakly 
directed bridges in the horizontal model. This algorithm is a Boltz- 
mann sampler That is, it involves a parameter x, and outputs 
a walk w with probability 

r <»» - WY 

where C{x) is the generating function of the class of walks under 
consideration. Of course, x has to be smaller than the radius of 
convergence of C. The average length of the output walk is 

The parameter x is chosen according to the desired output length. 
Boltzmann samplers have convenient properties. For instance, given 
Boltzmann samplers T_4 and Tg for two classes A and B, it is easy to 
derive Boltzmann samplers for the classes A + B (assuming ADB = 
0) and A x B. In the former case, one calls with probability 
A(x)/(A(x) + B(x)), and T B with probability B(x)/(A(x) + B(x)). 
In the latter case, the sampler is just (T^, Tg). If the base samplers 
run in linear time with respect to the size of the output, the new 
samplers also run in linear time. 

Moreover, if B C A, and we have a Boltzmann sampler for A, then 
a rejection scheme provides a Boltzmann sampler for B: we keep 
drawing elements of A until we find an element of B. 
Finally, if A = B x C, then we obtain a Boltzmann sampler for the 
class B by sampling a pair (b, c) and discarding c. 




FIGURE 14. Right: A random weakly directed bridge drawn with 
our algorithm. Above: zoom on a portion of the bridge. 



By Proposition [21 a weakly directed bridge is a sequence of partially directed irreducible 
bridges. We build our algorithm in four steps, in which we sample objects of increasing com- 
plexity. 

Step 1: The first step is to sample partially directed excursions. Let £ be the language 
of nonempty NES-excursions. As shown by Fig. 1121 this language is determined by the 
following unambiguous grammar: 

8 = E(l + E) + N£S + N£SE(1 + £). 
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We use this grammar to derive, first the generating function E(x) of partially directed 
excursions: 




, , 1 - x - x 2 - a; 3 - J(l- x 4 )(l-2x- x 2 ) 
E{X) = 2^ ' 

and then a (recursive) Boltzmann sampler for these excursions (see Section 3]). 
Step 2: The next step is to sample positive NES- walks, defined in Section 15. 1 1 More 
precisely, let Vn be the language of positive NES-walks that end with a N stepQ. We 
decompose these walks in the same way as bridges in Fig. [TSJ by cutting them after their 
last visit at height 0. We have the following unambiguous grammar: 

Vh = (N + EN+£EN)(1 +V N ). 

Given the Boltzmann sampler constructed for excursions in Step 1, we thus obtain a 
Boltzmann sampler for these positive walks. Their generating function is 

Pn{x) = 

Step 3: The object of this step is to sample irreducible NES-bridges; this is less routine 
than the two previous steps, as we do not have a grammar for these walks. To do this, 
we decompose the walks of Vh into irreducible factors (see Definition [14]) . Let Ze be the 
language of irreducible NES-bridges, and 7Z the language of irreducible positive NES- 
walks ending with a N step that are not bridges. Performing the decomposition and 
checking whether the first factor is a bridge or not, we find: 

V N =TZ + 2 E (1+V N ). 

We use this to construct a Boltzmann sampler for irreducible NES-bridges: Using a 
rejection scheme, we first derive from the Boltzmann sampler of Vn a Boltzmann sampler 
for the language Ie(1 + Vn)', these walks factor into an irreducible bridge, followed by 
a positive walk. We then simply discard the latter walk, keeping only the irreducible 
bridge. 

We construct symmetrically a Boltzmann sampler for the language Xw of irreducible 
NSW-bridges. We use another rejection scheme to sample elements of Zyv \ N. 
Step 4: Finally, the language VV of weakly directed bridges satisfies, as explained in the 
proof of Proposition [TBI 

W = l+I E W + (2w\N)>V. 
From this, we obtain a Boltzmann sampler for weakly directed bridges. 

Proposition 24. Let e be a fixed positive real number. The random generator described above, 
with the parameter x chosen such that xW (x)/W(x) = n, outputs a weakly directed bridge with 
a length between (1 — e)n and (1 + e)n in average time O(n). 

Proof. Let x > be smaller than the radius of convergence p of the generating function W, 
given by Proposition [THJ We first prove that if our algorithm outputs a walk of length m, it has, 
on average, run in time O(m), independently of the parameter x. 

The radius of convergence of the generating function Pn{x) is y/2 — 1, and is therefore larger 
than p. Hence the average length of a positive walk drawn according to the Boltzmann distri- 
bution of parameter x, being xP N (x) /Pn (x) , is bounded from above by pP N (p) / 'Pn(p) (see [TTJ 
Prop. 2.1]), which is independent of x. In particular, the algorithm described in Step 2 runs in 
average constant time (and the average length of the output walk is bounded). 



^One could just as well work with general positive walks, but it can be seen that this restriction improves the 
complexity by a constant factor. 
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Testing whether a positive walk is in 1Z can be done in linear time. Moreover, the probability 
of success in Step 3 is 

h(x) (1 + 1/Pn(x)) > x (1 + l/Pn{x)) = - (^(1- x i ){l-2x- x 2 ) + 1 - x + x 2 + .t 3 ) > y/2-1. 

Thus the average number of trials necessary to draw a walk of Ie (1 +Pn) is bounded by a constant 
independent of x. Therefore, the algorithm that outputs walks of Ie also runs in average constant 
time. The probability to draw in this step a walk of EN(1 + Vn) is x 2 (1 + 1/Pn(x)), which is 
bounded from below by x(\/2 — 1). Since in practice x will be away from 0, the probability to 
obtain an element of Ie distinct from N is bounded from below by a positive constant, so that 
we generate a walk of Ie \ N (or, symmetrically, of Iw \ N) in average constant time. 

Finally, the number of irreducible bridges in a weakly directed bridge of length to is less than 
to, so that the final algorithm runs in average time O(to). 

We now fix n and e, and choose x as described above (this is possible since xW'{x)/W(x) — > oo 
as x — > p). We call our sampler of bridges until the length m of the output bridge is in the 
required interval. Theorem 6.3 in [TT] implies that, asymptotically in n, a bounded number of 
trials will suffice. Indeed, the series W(t) is analytic in a A-domain, with a singular exponent — f 
(see Proposition [TS)) . I 

Fig. 14 shows a weakly directed bridge sampled using our algorithm, with a zoom on a portion 
of it. 
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